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Abstract 

The gyromagnetic relation - i.e. the proportionahty between the angular momentum L (defined 
by an inertial tensor) and the magnetization M - is evidence of the intimate connections between 
the magnetic properties and the inertial properties of ferromagnetic bodies. However, inertia is 
absent from the dynamics of a magnetic dipole (the Landau-Lifshitz equation, the Gilbert equation 
and the Bloch equation contain only the first derivative of the magnetization with respect to time) . 
In order to investigate this paradoxical situation, the lagrangian approach (proposed originally 
by T. H. Gilbert) is revisited keeping an arbitrary nonzero inertial tensor. A dynamic equation 
generalized to the inertial regime is obtained. It is shown how both the usual gyromagnetic relation 
and the well-known Landau-Lifshitz-Gilbert equation are recovered at the kinetic limit, i.e. for time 
scales above the relaxation time r of the angular momentum. 

PACS numbers: 
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The analogy between the dynamics of the magnetization in a magnetic filed on one hand, 
and the dynamics of a symmetrical spinning top in a gravitational field on the other hand, is 
often exploited in introductory courses on magnetism. The precession effect (i.e. the rotation 
of the extremity of a vector of constant modulus) is indeed easy to observe on a spinning 
top, while it is difficult to see with a ferromagnet because it would require observations at 
sub-nanosecond time scales. [l| However, the analogy seems to be incomplete because the 
dynamics of the symmetric spinning top implies inertial effects (e.g. nutation) while for 
a uniformly magnetized body, the dynamics of the magnetization is described by the time 
variation of the magnetization dM/dt (i.e. the velocity) and does not include the second 
derivative d'^M/dt^ (i.e. the acceleration). ^ In other terms, there is no inertia in the 
dynamic equation. The aim of this paper is to push the analogy to its logical end with the 
introduction of inertia [sl in the dynamics of uniform magnetization within the Lagrangian 
formalism. 

The precession of a uniform magnetic moment M = MgC^ {Mg is the magnetization at 
saturation and 63 the radial unit vector) under an effective magnetic field H is often presented 
as a consequence of the gyromagnetic relation M = 7L that links the magnetization to the 
angular momentum L. The constant 7 is the gyromagnetic ratio. The gyromagnetic relation 
and the value of the constant 7 = q/{2m) can be justified in a basic atomic model of an 
electron of charge q and mass m orbiting around a nucleus. This well-known model (see 
section IV below) constitutes the hypothesis of the Ampere molecular currents, validated 
by Einstein and de Haas in their famous experiments of 1915 - 1916. In the general 

case, with both spin and orbital contributions in condensed material, the gyromagnetic 
ratio writes 7 = gq/{2m) where the g factor accounts for the fact that the electron in a 
ferromagnet is a complex quasi particle, [l, 4| 

Using the gyromagnetic relation, the application of Newton's second law dL/dt = M x H 
leads directly to the precession equation dM/dt = 'jM x H. However the application of the 
Newton's law to a rigid rotating body (typically the spinning top in a gravitational field), 
leads to a more complex gyroscopic equation that contains inertial terms. As will be shown 
below, the gyromagnetic relation also imposes inertia for the dynamics of the magnetization. 
This paradoxical situation can be clarified by re-introducing the inertia in the equation of 
the magnetization and explicitly going to the kinetic limit. 

It is first useful to come back to the short history of the dynamic equation of the mag- 
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netization, especially with the introduction of the dissipation, since the precession equation 
dM/dt = 7M X H cannot account for the rapid relaxation toward the equilibrium state of 
the magnetization (typically after a couple of precession cycles, i.e. after some nanoseconds, 
in usual ferromagnets). In 1935 Landau and Lifshitz proposed an equation for the dynamics 
of the magnetization that takes into account both the precession and the relaxation along 
the magnetic field: dM/dt = 7M x H + h'M x (M x H), where h' is a damping term 
(defined below) and 7 = 7. I?! The basic argument used to derive the equation was to keep 
the modulus of the magnetization constant. The derivative dM/dt is hence perpendicular 
to the vector M. 

Two decades later, after the development of ferromagnetic resonance (FMR) experiments 
and motivated by the observation of systematic deviations from the above equation for 
high damping, T. L. Gilbert derived the equation that bears his name using a Lagrangian 
formalism. |9|, |lO| The dynamics of the magnetization is then described by the equation 



_=,Mx (1) 



with the introduction of the damping coefficient rj. 



The Landau - Lifshitz equation and the Gilbert equation are equivalent Uj provided 
that h' = and 7 = where a = jriMs is the dimensionless Gilbert damping (note 
that 7 7^ 7: this was the decisive improvement brought by Gilbert to the Landau-Lifshitz 
proposition). 

n 

In line with previous works performed by W. Doring, [2[ Gilbert introduced the La- 
grangian of a uniform ferromagnet with a kinetic energy T = LL : /~^/2, where / is the 
inertial tensor. He then chose an ad-hoc tensor of inertia in such a way that the inertial 
terms disappear from the dynamic equation (i.e. such that the Landau-Lifshitz equation is 
recovered at the low damping limit). To do that, a sufficient condition is to set to zero the 
two first principal moments of inertia Ii = I2 = (but keeping a non-zero kinetic energy: 



I3 7^ 0). As pointed out by Gilbert himself jlO| this puzzling condition does not seem to 
correspond to any realistic mechanical system (see footnote 7 : " / was unable to conceive 
of a physical object with an inertial tensor of this kind^). In the subsequent report about 
Gilbert's derivation, the ad- hoc and puzzling condition Ji = /2 = is explicitly stated 
despite its problematic character. In his presentation of the Gilbert equation published in 
1960 in the American Journal of Physics, ^] Brown wrote " We treat the rotating moment 



system as a symmetric top, with principal moments of inertia A = B = 0, C > 0. For a 
top made of classical mass particles, A = B = implies C = 0; but this top is not made of 
classical mass particles." In our notation A = Ii, B = I2, C = I3. In the reference textbook 
of Morrish fl3| (edited from 1965 to 2002), we can read: "A Lagrangian function, C, consis- 
tent with the accepted equation of motion (equation (10-3.2)) can he obtained by considering 
the magnetic system as a classical top with principal moments of inertia (0, 0, C)...". In our 
notation C = I3, and the equation (10-3.2) is dM/dt = 7M x H -\- damping. Accordingly, 
the mechanical approach is not presented as a realistic physical model (as it should, accord- 
ing to the gyromagnetic relation), but seems to be introduced as a pedagogical analogy of 
an unspecified non-classical theory, that would give a physical interpretation to the puzzling 
Gilbert's condition. Indeed, this strange condition is presented as a specific property of the 
magnetic moments that would be due to the fact that " this top is not made of classical 
mass particles" . 

This is probably the reason why, after more than half a century of intensive use of the 
Gilbert's equation, the full derivation following Gilbert's approach - with the complete set 
of principal moments of inertia (i.e. without ad-hoc assumption) - has not been proposed 

n n 

(see e.g. [ll|, |1J] for recent presentations of the Gilbert's derivation). However, as will be 
shown below, this derivation can be performed at an elementary level, as a direct application 
of the Lagrangian formalism. Although straightforward, this derivation is very instructive 
because it shows that the puzzling condition Ji = J2 = is not necessary to obtain the 
Gilbert equation. Instead, the Gilbert's condition is replaced by the necessary physical con- 
dition under which a diffusion process can be described by a non-inertial diffusion equation. 
This condition is the usual kinetic limit that results in the requirement that the typical 
measurement times should be longer than the relaxation time r of the momentum (here for 
the angular momentum r = Ii/ (qM"^)). 15| In this picture, the precession with damping is 
simply a diffusion process in a field of force, for which the angular momentum has reached its 
equilibrium. This change of paradigm has two consequences. A first important consequence 
is that an inertial regime of uniform magnetic dipoles is expected, and should be observed 
at short enough time scales. Second, the classical mechanical approach is much more than a 
pedagogical analogy, and it could be used (beyond the gyromagnetic relation) for a deeper 
understanding of non-equilibrium magnetomechanics and related processes. 
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I. THE MECHANICAL ANALOGY 



The mechanical model is sketched in Figure 1. A rigid stick of length Ms with one 
extremity fixed at the origin is described by the angles 9 and </?. The stick is precessing 
around the vertical axis at the angular velocity (p and is spinning around its symmetry axis 
at the angular velocity ij). The phase space of this rigid rotator is defined by the angles 
{9, ip, ■?/'} plus the angular momentum L. The relation between the angular momentum and 
the angular velocity Vt is L — IVt where / is the inertial tensor. 



A. The rotating frame 



In the rotating frame, or body-fixed frame {ei, 62, 63}, the inertial tensor is reduced to the 
principal moments of inertia {/i,/2,/3}. The symmetry of revolution imposes furthermore 

that Ii = I2: 



^ h ^ 



/i 
00/3 

In the fixed body frame, the angular velocity reads (see Fig. 1): 



v 



(2) 



Qi — (f sin9 sinip + 9 cosip 

0,2 — sin9 cosijj — 9 sintp (3) 
O3 — (p cos9 + ip 

— * 

The kinetic equation is obtained from the angular velocity: for any vector M of constant 
modulus carried with the rotating body, we have 



dM 



flxM 



(4) 



This equation can be inverted by cross multiplication by M and developing the double 

— * 

cross product. Since M — MgC^ we have: 



M dM ^ ^ 

= ^TTTT X — r— + ^363 



dt 



(5) 
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FIG. 1: Illustration of the magnetomechanical analogy of a spinning stick that precesses around 
the z axis. The coordinates of the stick in the space- fixed frame are parametrized by the angles 
{6,ip,ip) and the radius of the sphere is given by Mg. The body-fixed frame - denoted {61,62,63} 
- is spinning with the angular velocity ip and is precessing around Cz with the angular velocity ip. 

B. The Lagrange equation 

Following Gilbert and Doring, we introduce the Lagrangian of the system: 

where V{6, (p) is the ferromagnetic potential energy that defines the effective magnetic field 
H = — VV^. The effective field H comprises the applied field, the anisotropy field, the dipolar 
field (or the demagnetizing field), the magneto-elastic contributions, etc. 
The Lagrange equations are defined by : 

dtdqi dqi dqi 

The Qi refers to the three coordinates {6', yj, ?/^}, and the components of the kinetic mo- 
mentum are defined by the three derivatives = Lj. The function T is the Rayleigh 
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dissipative function. In a viscous environment, the Rayleigh function is defined by the 
damping coefficent r] such that T = ^{^^ = ^r/Mf (fi^ + 

For the magnetomechanical model, the Lagrange equations read: 



■j: [Ii9] — Iiip^ sinO cos9 + I^^tp sinO ( (f cos9 + ip 



av 

ae 



A 

dt 



Iiipsin^O + /3 ( cos9 + %l> ) cos9 



A 

dt 



I3 ( if cos9 + if) 



dip dip 
d4> 







(7) 



The right hand side of Eqs. ([7]) the last equation is equal to zero because there is no 
damping for spinning in the case of usual viscous environment. 1^ The quantity L3 = IsVls 
is then a constant of motion, and it can be written L3 = M^/'j without loss of generality. 



II. KINETIC EQUATION AND GILBERT'S ASSUMPTIONS 

It is not trivial to see how to recover the Landau-Lifshitz equation from Eqs. ([7]), even 
at the low damping limit 77—7-0. But it is clear that the inertial terms in the left hand side 
of Eqs. ([7]) are not welcome from that point of view and should be removed. The best way 
to consider the inertial terms is to take the kinetic Equation Eq. ([5]) with L = Ifl: 



MM dt ' ^ ^ 



(8) 



It is then rather immediate to see that the gyromagnetic relation L = M/'y cannot be 
recovered without removing the first term on the right hand side. This was the great idea of 
Gilbert to assume that Ji = 0. This is indeed a sufficient condition to kill the inertial terms, 
and the gyromagnetic relation is necessarily recovered with the definition 7 = Mg/L^ of the 
gyromagnetic ratio. 

With both assumptions, the Lagrange equations Eqs. ([7]) rewrite: 







A 

dt 



Ms 



cos6 



dip 



9V_ 

dip 



(9) 



A 

dt 



Ms 
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Since in the rotating frame the effective field H = — VV^ reads {Hi = H2 = 

%,}, inserting Eq. @ into Eq. leads to : 

— ^,M.[H-,—j (10) 

This is the well-known Gilbert's equation Eq. ([1]) obtained following the standard La- 
grangian approach. 



However, the absence of inertia shows that the equation should be derived in the config- 

n 

uration space instead of the phase space (this is performed e.g. in references p/Z|)- Indeed, 
the dynamics is described by the two variables 6 and (f and not in the phase space defined 
by the five variables 6, ip and the components of L. Accordingly, the gyromagnetic relation 
is - in this approach- not necessary (nor sufficient) for the derivation of the Gilbert equation. 



III. BEYOND GILBERT'S ASSUMPTION 



In order to take into account the gyromagnetic relation, it is necessary to go beyond the 
Gilbert's ad- hoc assumption and to set Ii = I2 ^ 0. The generalization of the Ampere 
molecular model from a quasi-one dimensional atomic model (the electric charge q of mass 
m distributed along the circular orbit) to a more realistic three dimensional atomic model 
- for which the orbits form an ellipsoid of revolution (the electric charge is now distributed 
in three dimensions) - imposes non- vanishing inertial moments Ii = I2 Indeed, 
two parameters are necessary to take into account the amplitude of the magnetic moment 
(Ms = 7/(^3/3)) on one hand, and the anisotropy of the ferromagnetic material (with the 
dimensionless parameter 1 — /1//3) on the other hand. 

Note that for the magnetic system, the variables ip and ip are not 
removed from the model. Let us take for the sake of simplicity iJj = 0. 
L3 = ^, Eq. dZD gives: 

(l2 = ^3 (1 - ^) ^1 + fHi (11) 

where the typical relaxation time r = has been introduced. The relaxation time r 



efined and should be 



18 1 With the relation 
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is the typical time above which the diffusion approximation is vahd, i.e. above which the 
angular momentum has relaxed toward the equilibrium state, jsl 

Using Q.e^ = 0, Eq. ([1]), Eq. (E]) and the time derivative of Eq. (E]), Eq. ( ITT]) re- writes: 



dM 



jM X 



, dM d'^M 
H -r]{ — + T- 



(12) 



dt ' ' dt"^ 

This is the Gilbert's equation of the dynamics of the magnetization that includes the new 
inertial term — 7rM x [d^M /dt^ 



IV. TYPICAL TIME SCALES 

The limit of Eq. (1121) for t r, where r = leads to the kinetic limit. Since the 
damping 77 can be replaced by the usual dimensionless Gilbert coefficient a = 'jrjMs, we have 
T = tJ-^^- a rigorous study of the asymptotic behavior (as a function of the parameters 
'jH, and 77) is beyond the scope of this work. However it is sufficient to observe that 
the limit r — > leads straightforwardly to the LLG equation: 

dM - ( -^^ dM\ 

In the same manner the vectorial gyromagnetic relation is recovered at the limit r — > 0. 
Eq. (IHI) gives: 

dL ( ^ d^M\ IdM IdM 

— = r]T M X — — + > — (14 

dt ' \ dt^ j dt dt ^ ^ 

The sufficient condition of validity of the Gilbert equation Ji = is hence replaced by 
the condition r — )■ (i.e. t <^t for the relevant range of the parameters). 

An estimation of the value of r gives the typical time scale for which inertial effects can be 
observed. Here we come back to the simplest argument for the justification of the value of 7, 
namely the model of the Ampere molecular currents. This is a quasi-one dimensional atomic 
model, for which the atomic orbital moment is defined by the electronic charge q orbiting 
around a nucleus at the distance r with a velocity v. This system defines an electric loop 
that generates a magnetic moment M = IS Cz, where X = qv/{27ir) is the electric current, 
S = vrr^ is the surface enclosed by the loop, and Cz is the vector normal to the loop. This 
leads to the microscopic magnetic moment Mg = qvr/2. If we take the Bohr radius Qq and the 



electron velocity v with the Heisenberg relation mvao > h/2, the Bohr magneton is obtained 
for the minimum value of the atomic magnetic moment ij,b — 7^/2. On the other hand, 
the angular moment of this system is L3 = rmv and the ratio M^/L^ = 7 = q/i^m). The 
angular frequency Q3 is given by L3 = /a 1^3 = fJ,B/l, i-e. ^3 = fis/i'lh) where I3 = mal. 
We have ^ 3 10^^ rad/s and an order of magnitude of the typical times r = ^'^^ 
which inertial effects should be observed is around a femtosecond (for a damping coefficient 
a such that Ii/{l3a) ^ 0.1). 

V. CONCLUSION 

The paradoxical role played by the angular momentum for the dynamics of the mag- 
netization has been studied in the light of the model introduced by T. H. Gilbert for the 
demonstration of the equation that bears his name. The demonstration has been recon- 
sidered without the puzzling Gilbert's assumption of vanishing first moments of inertia 
Ii — I2 — and Is 7^ 0. Instead, a general inertial tensor with the three arbitrary principal 
moments of inertia {7i,7i,73} has been used. A generalized expression of the equation of 
the dynamics of the magnetization is obtained, that includes an inertial term: the mechan- 
ical analogy of the magnetic moment with the rigid rotator is complete. Both the usual 
expression of the Landau-Lifshitz-Gilbert equation and the usual gyromagnetic relation are 
recovered provided that a kinetic limit is performed for time scales much larger than the 
relaxation time of the angular momentum r = Ii/{r]M'^). The typical time scale is found to 
be of the order of the femtosecond. 
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